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Abstract. Experimentally, the production of oppositely charged dimuon events by neutrino and anti- 
neutrino deep inelastic scattering (DIS) is used to determine the strangeness asymmetry inside a nucleon. 
Here we point out that the direct production of D-meson in DIS may make substantial influence to the 
measurement of nucleon strange distributions. The direct D-meson production is via the heavy quark re- 
combination (HQR) and via the light quark fragmentation from perturbative QCD (LQF-P). To see the 
influence precisely, we compute the direct D-meson productions via HQR and LQF-P quantitatively and 
estimate their corrections to the analysis of the strangeness asymmetry. The results show that HQR has 
stronger effect than LQF-P does, and the former may influence the experimental determination of the 
nucleon strangeness asymmetry. 



1 Introduction 

Studying strange and anti-strange quark distributions of 
a nucleon is an important part in the study of the nucleon 
structure. An asymmetric strange distribution, i.e., the 
parton distribution function (PDF) of strange quark being 
not equal to that of anti-strange quark inside a nucleon, is 
naturally predicted by some non-perturbative models [H 
ElE]. Further clear check of the strangeness asymmetry is 
not only important for the study of the nucleon structure 
itself, but also for understanding relevant phenomenon 
in some experiments. For example, the so-called NuTeV 
anomaly phenomenon [?J[5] can be explained by nonzero 
strangeness asymmetry pl[7Il8ll Wll)l[Tl"] . 

Because of the smallness of the strange and anti-strange 
components in a nucleon, to measure the strangeness asym- 
metry is a challenging job indeed. The most sensitive re- 
action to measure the strange and anti-strange distribu- 
tions is the production of dimuon events in neutrino and 
anti-neutrino nucleon deep inelastic scattering (DIS). To 
leading order (LO) of the dimuon production, the events 
are caused by the charged-current (CC) charm production 
subprocesses + s (d) — > yT + c or + s (d) — > /i + + c 
and a cascade decay c-^y + + -- -orc-> u~ + ■ ■ ■. The rel- 



evant transition v^v^) 



[i + c or v 
-> /i + + • • • or c - 

+ d(d) — ► n~ (/i + )+c(c) is Cabibbo 
suppressed whereas ^(77^) + s(s) — > ^(p + ) + c(c) is 
Cabibbo favored. Thus the dimuon events are sensitive to 
the strange and anti-strange distributions of the target nu- 
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cleon. In the literature, measurements of the strangeness 
asymmetry via oppositely signed dimuon are carried out 
by CCFR and NuTeV [T2 l fT ^ T4 T fl3 1 [l ^ [T7] experiments. 

The results of global analysis [TSjEllTS] indicate the 
strangeness asymmetry S~ = J £[s(£) — likely to 

be positive, e.g., in Ref.pU -0.001 < S~ < 0.005 is ob- 
tained. The recent NuTeV reanalysis up-to next-to-leading 
order (NLO) of perturbative QCD (pQCD) with improved 
method supports the positive strangeness asymmetry S~ = 
0.00196 ± 0.00046(stat) ± 0.00045(syst)±8;$£$ 
(external) [T7]> which is consistent with the global anal- 
ysis, although at early stage the analysis of CCFR and 
NuTeV dimuon events at LO and even NLO do not sup- 
port the strangeness asymmetry [T3l[T4l[T5] . 

In this work, we take a systematic study on the in- 
fluence of direct D-meson production at order a 2 s to the 
determination of the nucleon strangeness asymmetry via 
dimuon events in neutrino experiments. With consider- 
ation of the experimental kinematic cuts in CCFR and 
NuTeV, we point out that there are two kinds of direct 
D-meson production: heavy quark recombination process 
(HQR) and light quark fragmentation in the pQCD pic- 
ture (LQF-P), which can contribute to the cross section 
difference between neutrino and anti-neutrino induced CC 
DIS. The direct D-meson production may influence the 
strangeness measurement (determination) that depends 
on their magnitude, so we further calculate the produc- 
tion quantitatively and investigate their influence to the 
strangeness determination. Although our preliminary re- 
sult on HQR process has been briefly reported in Ref . [20] , 
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in this paper, we would present it in more detail and with 
some improvements. From final results we find that the 
influence of the direct _D-meson production to the mea- 
surement of the strangeness asymmetry could not be neg- 
ligible. 

The rest of the paper is organized as follows. In sec- 
tion II, we discuss CC charm production to the dimuon 
cross sections at LO and NLO. And then we present the 
two kinds of direct Z?-meson production: HQR and LQF- 
P, which in fact are of high order ones, and show how 
the two processes can affect the extraction of the nucleon 
strangeness asymmetry. In section III, we show the numer- 
ical calculation about HQR direct production of D-meson 
and estimate the influence due to the direct production. In 
section IV, we show numerical calculation and estimate as 
in previous section about those for LQF-P production. In 
conclusion section (section V), we summarize and discuss 
the obtained results. 



the differential cross section can be expressed as 



d 2 a 



v^A—tfi-fi+X 

d£dQ 2 



r 2 

- F t 

~ 9 J< 



Pd(tQ 2 ) + Nu(d,Q 2 )^ |2 

\ V cd\ 



P + N 



-s(Z,Q 2 )\v cs f 



D, 



ef 



(2) 



where the PDFs in the neutron is related to PDFs of pro- 
ton by dn(g) = u(f), s„(£) = s(£) etc. 

Since \V cd \ 2 ~ 0.05 and |y cs | 2 - 0.9 [22], the v„ in- 
duced dimuon cross section is sensitive to the strange dis- 
tribution s(£, Q 2 ) in the target nucleus. Similarly, in- 
duced dimuon cross section is sensitive to anti-strange dis- 
tribution s(£, Q 2 ). The difference between dimuon cross 
sections induced by and F M is directly related to the 
strange distribution asymmetry. The difference to LO can 
be expressed as: 



d 2 a 



2 Dimuon events and direct D-meson 
production in neutrino DIS 

Experimentally with the CCFR and NuTeV detectors, the 
oppositely charged dimuon signal induced by CC charm 
production in (V^) DIS has a distinct feature and is 
not very difficult to be detected. The first muon of the 
dimuon is from the v^{p^) vertex, and the second muon is 
from a little delayed muonic decay of the produced charm. 
The life time of it and K is much longer than charmed, 
so those muons from tt or K meson decay can be largely 
eliminated, i.e., they will not contribute to the dimuon 
events concerned here. 

According to pQCD factorization theorem, for ^-proton 
DIS, to LO the differential cross section for dimuon pro- 
duction induced by CC production of charmed hadron H 
can be expressed as [21] : 



d£dQ 2 



d 2 o- 



d£dQ 2 



-JLf c B ef {[ s (tQ 2 )-nt,Q 2 )]\V cs \ 2 

+ pTn t Pdvit q2) + Nu *& q2) ] |ycd|2 } ' (3) 



d 3 c 



i+X 



G 



d£dQ 2 dz 



^fc[d(^Q 2 )\V cd \ 2 + s(tQ 2 )\V cs \ 2 ] 

w 

Y,D?{z)Br H . (1) 



7TT 
X 



H 



where Q 2 ) and s(£, Q 2 ) are the parton distribution 
functions (PDFs) of d and s quarks in the proton, and 
£, relating to the Bjorken scaling variable x through £ = 
x(\ + m 2 /Q 2 ), is the light-cone momentum fraction of the 
struck quark; Q 2 = — (p — k) 2 is the minus squared invari- 
ant momentum transfer with p and k being the momen- 
tum of the incident and the scattered pT respectively. 
r w = 1 + Q 2 /M^ and f c = 1 - m 2 c /S(, with M w being 
the W-boson mass and S being the squared CM. energry 
of the neutrio- proton system; Df? (z) is the fragmentation 
function for a charm quark to the charmed hadron H, and 
Brn is the branching ratio of muonic decay for H . Carry- 
ing out the integration over z and the summation over H, 
with the definition B e f = J dz^ H {z)Brn 1 for target 
nucleus with proton number P and neutron number N, 



where d v (^Q 2 ) = d(t,Q 2 ) - d(£,Q 2 ) and u v (^,Q 2 ) = 
M (£> Q 2 ) — Q 2 ) are valence distributions of proton. 

The Feynman diagrams for CC charm production at 
NLO from subprocess oiv^N DIS are shown in FIG. 1 |16j . 
In fact, the last two diagrams of FIG. 1 [16] are involved in 
the leading logarithm (LL) evolution for the parton flavor- 
singlet components of PDFs, and the flavor singlet compo- 
nents of the PDFs contribute to the dimuon cross sections 
symmetrically for the v^- and F^-induced DIS. Thus even 
up-to NLO, the difference between v^- and j/ M -induced 
dimuon cross section is still proportional to [s(£, Q 2 ) — 
s(£,g 2 )]|^ s | 2 +p^v[P^(£,Q 2 ) + iV^(£,Q 2 )]|K d | 2 ,i.e., 
we need consider neither the last two diagrams of FIG.l 
nor the flavor singlet components of PDFs in the dimuon 
cross section difference. 

However, the direct ZJ-mesorQ production that convo- 
lutcs to the valence components in the nucleon, as we will 
discuss in the following, can raise the rate of dimuon pro- 
duction after experimental kinematic cuts, therefore de- 
terminations of strangeness asymmetry can be distorted 
by the direct £>-meson production in certain degree. 

Now let us focus the light on the contributions from the 
direct D-meson production, although according to pQCD 
the lowest order Feynman diagrams of the direct produc- 
tion are of the order a 2 . 



The D* meson, the excited (cq) bound state in 3 Si, has 
very great cross section in production to compare with the D- 
meson production, and decays into D meson via strong and/or 
electromagnetic interaction with almost 100% branching ra- 
tio, thus the consequence of the D*-meson production will be 
as direct D-meson production in accounting the muons in the 
dimuon events. Therefore throughout the paper "D-meson pro- 
duction" always mean the production of D and D* . 
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(b) 




Fig. 1. Feynman diagrams for calculating the charm produc- 
tion by energetic collision via CC weak interaction up-to 
NLO. The thick lines denote those of c or c quarks. 





(c) 



(d) 



One of the direct £)-meson production mechanism, the 
so-called heavy quark recombination (HQR) process, at 
LO level is described by the diagrams in FIG. 2. It is 
stimulated from the heavy quark recombination mecha- 
nism 23, 24, 25J , which combines a heavy quark and a light 
anti-quark of similar velocity to form a meson. Refs. [24l 
125) employ simple pQCD pictures and explain the charm 
photoproduction asymmetry and the leading particle ef- 
fect [26] successfully. 

Namely the difference between v^- induced and v^- 
induccd dimuon cross sections caused by the HQR pro- 
duction of direct D-meson can be computed: 



didQ 2 



didQ 2 



HR 



= J2j dx[q(x,^)-q(x,^)}^^Br D , (4) 

where q denotes a possible light quark in the target and D 
denotes a produced £>-meson. dajy denotes the differential 
cross section of the subprocess 



+ q -> fj, +s(d) + D, 



(5) 



which, in terms of CP transformation, is known to be 
equal to that of the subprocess 



+ q -> M + + s (d) + D . 



(6) 



And Btd denotes the muonic decay rate of the D-meson. 
In fact, from Eq.(j4]), it is easy to realize that only the 
u — d valance components of the PDFs contribute to the 
cross section difference (s — s contribution comparatively 
negligible), thus q could be either u or d. Such a con- 
tribution may reduce the strange distribution asymmetry 
S-(£,Q 2 ) ee £[s(£,Q 2 ) - s(£,Q 2 )] value obtained by the 



Fig. 2. (a) and (b) are diagrams for (cq) recombination into a 
D meson in induced process Eq. ([5j ; (c) and (d) are diagrams 
for (cq) recombination into a D meson in 77 M induced process 
Eq.©. The thick lines in the diagrams denote those of heavy 
quarks, and the shaded blobs denote D or D mesons. 



analysis where the HQR process has been ignored. Thus, 
roughly speaking, the true value of the strange distribu- 
tion inside a nucleon S~ '(£, Q 2 ) may be changed from the 
existent analysis result which has ignored the HQR effects 
by a positive correction factor SS^ K (^, Q 2 ) 



S^t, Q 2 ) = Sl naly & Q 2 ) + SS^fo Q 2 ) 



(7) 



For a quantitative estimate of the HQR correction to 
the measured strangeness asymmetry, we suppose that the 
Cabibbo suppressed valence contribution in Eq.® has 
been deducted in LO (or NLO) analysis of dimuon events 
though present knowledge of the PDFs, 



s(f Q 2 n LO = \ - dl ° V » 

, Pd v (t;,Q 2 )+Nu v (Z,Q 2 ) 

P + N 



ho\V cs \ 2 [s(C,Q 2 



-fho\V c d\ 



(8) 



Where / L o = G|./ c i?r c /7rr 2 , is the coefficient in Eq.Q 
for LO cross section, and 



d£dQ 2 



is the differential cross 



section for ^(F^-induced dimuon production. While in 
fact, the measured cross section also includes the contri- 
bution from HQR as Eq. ^ , which should be deducted to 
obtain the real strange distribution asymmetry 



ho\V cs \ 2 [s(^Q 2 )-s(tQ 2 )} 



real 



d 2 G Vu 



didQ 2 d£dQ 2 



d 2 <Jv» 

d£_dQ 2 



d£dQ 2 



HQR 
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Fig. 3. Diagrams for the LQF-P process in induced CC DIS. 
The produced D meson will decay partially into fi + to form 
the second muon. Thick lines denote that of heavy quarks, and 
shaded blobs denote the D meson produced directly. 



ho\V cd \ 2 [Pd v (^Q 2 )+Nn v (^Q 2 )] 



P + N 



fix)\Vc,\WZ,Q i )-'3(t>Q*)]£Ly 



d£_dQ 2 d£dQ 2 



(9) 



HQR 



LQF-P production can be expressed as 



d 2<J v u N^u,- u+X 



d£dQ 2 



LQF-P 



E 



dx Pu(x)+Nd( X ) da-^ BrD 



P + N 



d£_dQ 2 



Where d&cD denotes the cross section of subprocess + 
u — > + D + c (diagrams in FIG. 3). The F M -induced 
dimuon cross section from LQF-P production can be ob- 
tained by CP transformation to that of v, and thus the 
Vfi and T/^-induced dimuon cross section difference from 
LQF-P process can be expressed as 



N^u+u-X 



d£dQ 2 



d£dQ 2 



J LQF-P 



Pu v (x) +Nd v (x) ^ da-cD D ,..„, 
dx — 2 , Jt - J ^ Br n- ( 13 ) 



P + N 



d£dQ 2 



Thus, the real strange distribution asymmetry S (£, Q 2 ) 
f i Q 2 ) ~ Q 2 )] can be deduced. 



^-LO 
analy 



(Z,Q 2 ) 



r ^2 



/LolKsl 5 



cTer,, 



d 2 aiy.. 



d£dQ 2 d£dQ 2 



(10) 



HQR 



Thus for LO analysis the HQR correction to strange dis- 
tribution asymmetry can be estimated from Eq. (l4|7|10|) 



«5S HR (£,Q 2 )^ G 2 f ^b W J\ Vcs \2 

x E/ dx[q(x,/j, 2 ) - q(x,fi 2 

q,D 



, d 2 a D 
1 d£dQ 2 



Br D (11) 



For NLO analysis, the factor /lo m Eq. ([T0[) should be re- 
place by a NLO coefficient, and the HQR correction should 
differ from this LO one Eq. fTTj) by a K factor. Numerical 
calculation of <55 I I J R (^, Q 2 ) will be the content of the next 
section. 

The other possible mechanism is the so-called light 
quark fragmentation process in the order a 2 of pQCD 
(LQF-P). Its subprocess is described by FIG.3. From FIG.3 
it is easy to realize that LQF-P contributes not only to the 
inclusive oppositely charged dimuon events but also to 
trimuon and the inclusive same charged dimuon events. 
Here we focus only on the contribution from LQF-P to 
the inclusive oppositely charged dimuon with the restric- 
tion that the second muon comes from the produced direct 
Z?-meson. This mechanism is considered as fragmentation 
of a light quark into a D-meson in perturbative picture, 
which is different from the non-perturbative picture dis- 
cussed in Ref. [27] . 

For nucleus target with proton number P and neutron 
number N, the i^-induced dimuon cross section from the 



Note here that to compare with the valance components, 
the other components are tiny, so in Eq. (|13p . except the 
valance components, the other components are ignored. 

Numerical calculation of the LQF-P process will be 
presented in section IV. 



3 Calculation of the HQR process 

As pointed out in the above section, the direct D me- 
son production from HQR process (diagrams in FIG. 2) 
may influence the measurement of the nucleon strangeness 
asymmetry. Our calculation of the HQR process follows 
the method in Rcfs. [23 , 24, 25j . Namely we 'factorize' the 
production into the two steps, one of them (the one we 
start with) is the production of the relevant 1 S'o or 3 Si 
state of (cq), the other one (the following one) is the 
combined quark pair to evolve into the relevant D- or 
£)*-meson under certain possibility, that is determined by 
other experiments. 

When a charm c quark and a q light anti-quark with 
momentum p c and p q respectively are produced with p q 
to be small in the c-quark rest frame, the c and q can be 
constructed into a So color-singlet state (cq) if q- and 
c-quark have the same color, then we have the following 
substitution in amplitude for the combination of the par- 
tons: 

6 

Vj(p q )ui(p c ) — > x q j^-m c {^ c - m c )7 5 , (14) 

with p q = x q p c and limit x q — * may be taken. Ac- 
cordingly, for the color-singlet 3 Si state (cq), 75-matrix 
in Eq. (fT4]l should be replaced by ft, where e fl (P,X) is 
the polarization vector of the 3 Si (cq) state with total 
momentum P and polarization A of the quark pair. 5ij 
in Eq. (fl~4")) is the color factor for color-singlet state with 
quark colors i ,j = 1,2, 3. Corresponding to the diagrams 
in FIG. 2(a) (b), the amplitude for the production of the 
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(cq) color-singlet 1 5'o state induced by may be written 
down directly: 



M 



16tt GFOL s m c 5 n 



9V2r w (2l- Pc ) 
x u(l)~f"(j> c - m c )75 



V cf u(k)j»(l-j 5 )u(p) 
j> - - ft s + m c 



( P -k-k s y 



■7m (1 



+7^(1-75) 



v(k s 



- 75) 

(15) 



Here p, I, k and k s denote the momenta of the , the 
initial incoming q, the produced \x~ and the produced s 
or d respectively. V c f is the CKM matrix element with 
/ = s, d. Color factor T^SijT^ — |<5 mn has been included 
in the amplitude of Eq . (fTS"]) . Squaring the amplitude, aver- 
aging over the spin and color of the particles in the initial 
state, and summing up the spin, color and possible flavors 
(/ = s, d) of the particles in the final state, the 'averaged 
and summed' squared amplitude is 



\M\ 2 [(cq)l So ] 



4ir 2 l6 3 G 2 F a 2 m 2 c 

~ 81r2 £( s _ A _ £)3 

x |(7 2 [(s - A)m 2 c + (s- A - B)Q 2 } 
+ m 2 c C[(s - A)m 2 c -(s-A- B)(s - 
+ m 2 c A[s(Q 2 + m 2 ) -BQ 2 ]}. 



n 

(16) 



Here s = (p + I) 2 = Sx is the squared CM. energy of 
the subprocess. The variables A, B and C are defined as 
A = 2(1- k), B = 2(1 ■ k s ), C = 2(p- k s ). 

For the color-singlet 3 5i (cg)-state production, the 'av- 
eraged and summed' squared amplitude is 



\M\ 2 [(cq)l Sl ] 



4ir 2 16 3 G 2 F a 2 m 2 c 



81rlB(s -A- BY 
x {3C 2 \(s - A)m 2 c + (s - A - B)Q 2 ] 

+ C[(3Q 2 - 5s + 2A)(s - A- B)m 2 c 

+ 3(s - A)m\ + A(B — S)(s — A — B)Q 2 } 

+ [s(Q 2 + m 2 ) - BQ 2 ] 

x [2(s -A- B)(s -B- m 2 c ) - Am 2 ] Vl7) 

For the color-octet (cq) either in 1 So or 3 Si state, the 
amplitude differs from the 'color-singlet' ones only in the 
color factor, e.g., the S mn in Eq. (fT5|) should be replaced 

by with a = 1,2,..., 8. 

Then the produced (cg)-states will evolve into either D 
or D* meson with certain probabilities. The cross section 
for the subprocess of a Z?-meson production + q — > 
pT + D + s (d) can be expressed as 



da D = ^2 da [cl!)% p((cq) c s -> D) 



(18) 



c,s,q 



Where c and s denote the color and angular momentum 
quantum numbers of the (cq) state, and p((cq)g — > D) is 



the non-perturbative parameter characterizing the proba- 
bility for the (cg)g to evolve into a state including the D 
meson. 

For (cq)g with various antiquark q, spin (s) and color 
(c), there may be a number of p parameters. However, 
the number can be greatly reduced in terms of symme- 
tries and some approximations. First of all, according to 
the color-factor of the processes, the cross sections for the 
production via color-octet and color-singlet (cq) differ by 
a single factor | , therefore we can express the cross section 
for D- meson production as follows: 



d(T D = ^2 d ^(cq)i PeS((cq)s -> D) 



(19) 



with the definition: 



Pcs((cq) s —* D) = p((cq)l —> D) + -p((cq)t D) . 

If the produced D-meson has different flavor from that of 
the quark pair (cq), that means the quark pair (cq) must 
emit a flavored object (such as a pion etc) in the mean- 
time forming the D meson, e.g., (cu) — > D + + ir~ , then 
the relevant p e g will be relatively suppressed in the large 
N c limit. As in Refs. [24ll25j . we neglect such transitions 
from the quark pair (cq) to a different flavored D-meson. 
Furthermore, SU(3) light quark flavor symmetry indicates 
Pctt((cu) s — ► D a ) ~ p c fi((cd) s — > D + ). As discussed in 
Ref. [24] . heavy quark spin symmetry implies 

p cS ((cq)i So - D(cq)) = p cS ((cq) 3Sl -» D* (cq)) . (20) 

Thus, only two independent parameters are left: 

Psm = Pes((cd)is -> D + ) = p cS ((cd)3 Sl -> D* + ) 

= PeS ((cu) 1So -> D°) = PeS ((cu)s Sl -> D*°), (21) 

Psf = Pctt({cd)i So -> D* + ) = p cB ((cd)3 Sl -> D + ) 

= p eS ((cu) 1So - D*°) = Pctf ((cu) 3Sl - D°). (22) 

Now 5Sft R (£, Q 2 ) caused by HQR, according to Eq. lfTll 
can then be evaluated by means of Eqs. (|16|17p with the 
auxiliary parameters p sm and p s f which are determined 
from relevant experiments. For the target to be proton, 
the HQR correction of Eq.lfTTj) can be expressed as 



SS HR (^Q 2 ) 



11K ' ' " G 2 F f c B ef \V cs \ 2 

, I <'x I [d v (x, p 2 )Br D + + u v (x, p 2 )Br D a] 



x 



d 2 b 



-Ps 



( C 9)3 S , 
2~P*f 



didQ 2 rs "" ' d£dQ 
+ [d v (x, p 2 )Br D ,+ + u v (x, p 2 )Br D ,(i] 

( d2& (cq)i So , d2a (cqh Sl \\ foo , 
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where the subprocess cross sections do"( C g) ls and da^ c ^ 3s 
are related to the 'averaged and summed' squared ampli- 
tudes of Eq. (|16I17[) and are independent of quark flavor 
q. And for the nucleus target with proton number P and 
neutron number N, the HQR correction can be expressed 
as 



G 2 F f c B ef \V C! 



Pd v (x,n 2 ) + Nu v (x,n 2 ) 



P + N 
Pu v (x,/j, 2 ) + Nd v (x,n 2 ) 



P + N 



( 



d 2 a 



+ 



d£dQ 2 r '" 1 ' d^dQ 
Pd v (x, y 2 ) + Nu v (x. [i 2 ) 



Br D o 
i P»f) 



P + N 
Pu v (x,ij, 2 ) + N d v (x,n 2 ) 



Br 



D* + 



P + N 

—Psf + 



Br D *a 



d£dQ 



d£dQ 2 



Ps 



0}(24) 



For a quantitative estimate, in the following, we cal- 
culate the <5Shr(£' Q 2 ) f° r isoscalar target (P = N), with 
the parameters p sm = 0.15, p s f = 0, which are taken from 
the extraction from experimental charm photoproduction 
asymmetry by Ref . |24) . In this case, the HQR correction 
of Eq. (j2"4|) can be simplified, 



G 2 F f c B ef \V C! 



dx 



d v (x, fj 2 ) + u v (x, /i 2 ) 



(Br D+ +Br D o)- 



d 2 a 



( c g)i 



-{Br D ,+ +Br D ,o)- 



d£dQ 2 
d 2 a, 



d£_dQ 2 



(25) 



In fact, for calculating Eq. ([25f wc need to carry out 
a three-dimensional integration numerically. In the CM. 
frame of the subprocess, £ and Q 2 are related to the energy 
ko and the angle 9\ of the outgoing muon (from or V ^ 
vertex) relative to the incident direction: 



i = 



x[y/sko(l — cos^i) + ni 2 
s - \fska{l + cos^i) 



(26) 



Q 2 = y/sko(l - cos6»i), 



(27) 

Bearing the NuTeV dimuon experiment in mind, the 
incident energy of neutrino or anti-neutrino is taken to 
be E v — 160 GeV, which is approximately averaged value 
of the experiment. Furthermore, for our numerical calcu- 
lation, charm mass m c is fixed with the value 1.5 GeV, 
the coupling constant a s {n) is running as specified in 
CTEQ6L[28] and the parton distribution functions of the 
nucleon are taken from CTEQ6L[28] too. The branching 
ratio for the muonic decay of D meson and the CKM ma- 
trix elements are taken to be the central values in Ref. [22] , 
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Fig. 4. SS HR as a function of £ for Q 2 = 10 GeV 2 (solid 
lines), Q 2 = 20 GeV 2 (dash-dotted lines), and Q 2 = 30 GeV 2 
(dashed lines) . Thick lines are results for /i = /io = \Zp 2 ± + rn^ 
and thin lines are for jj, — Q. 



and the B e f is taken to be the central value of B c _>„+x in 
Ref. 22J. Since the two opposite charged muons in NuTeV 
experiment are required to have energy greater than 5 
GeV, so we make a cut for the produced /i and the D 
meson accordingly. 

The obtained result of as a function of £ for 

Q 2 = 10 GeV 2 (solid lines), Q 2 = 20 GeV 2 (dash-dotted 
lines), and Q 2 = 30 GeV 2 (dashed lines) is shown in FIG. 4. 
Since the calculation is of LO for the direct D-meson 
production, there are theoretical uncertainties, such as 
that from the energy scale fj, of perturbative QCD, so 
to see the uncertainty, we calculate with two types 

of choices about fi. The thick lines in the figure present 
the results where the factorization scale /i is taken to be 
/i = \/ v\\_ + m c = A*0) where p c ± is the transverse mo- 
mentum of the produced D meson to the direction of the 
W boson in the nucleon rest frame. It is an analogous 
choice to that in Ref. [24] , where its relevant charm pho- 
toproduction is calculated under the factorization scale 
/i = \/pj_ + m 2 with p± being the transverse momentum 
of the produced D to the incident photon direction. The 
results in an alternative scale choice fj, — \J Q 2 — Q, are 
also shown by the thin lines in FIG. 4. One fact that should 
be noted here is that the cross section from the HQR pro- 
cess decreases very slowly with the increase of the energy 
cut of the produced D meson E cu t, namely, the 'recombi- 
nation' is not suppressed by the cut taken in experiments 
very much. That can be understood by the fact that the 
difference of the direct production of D-meson by the HQR 
is related to the subprocess with valance quark inside a nu- 
cleon (Eq . (fTT]) ) so that the D-meson relevant to the dif- 
ference can carry comparatively high energy (momentum) 
that escapes from the cut quite a lot. 

From FIG. 4 one can see that at fixed Q 2 in each case, 
<55jj R (£) peaks in the region £ = 0.1 — 0.2, over the peak 
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Fig. 5. as a function of Q 2 for £ = 0.06 (solid lines), 

£ = 0.15 (dash-dotted lines), and £ = 0.3 (dashed lines). Thick 
lines are results for /i = fio — yP%± + m l an d t nm lines are 
for [i = Q. 



Fig. 6. ^S'jjr as a function of m 3 for £ = 0.15 (solid lines) 
and £ = 0.3 (dashed lines) at Q 2 = 20GeV 2 . A plot to show 
the resultant uncertainties from choosing m s in treatment of 
the collinear singularity. Thick lines are results for /i = fio — 
\/p 2 ± + m% and thin lines are for fi = Q. 



<5Sg R decreases with Q 2 increases, and the results with 
factorization scale [i = Q are smaller than those with /i = 
/iQ. The uncertainty from the choice of the factorization 
scale [i can also be seen when the scale /i = /^o is varied 
by a factor of 2: the results become nearly trebles when 
/i = /io/2, and the results reduce nearly by half when 
/i = Generally the uncertainty may be suppressed 

by NLO calculation, but we leave the study beyond the 
present calculations. 

The behavior of SS^ R as functions of Q 2 is shown in 
FIG. 5: the solid lines are those for £ = 0.06, the dash- 
dotted lines are those for £ = 0.15 and dashed lines are 
those for £ = 0.3. The thick lines are results for factoriza- 
tion scale fj, = [1q and the thin lines are those for fi = Q. 

In our calculation, a colinear singularity may arise from 
the strange quark propagator in diagrams FIG. 2(b) and 
FIG. 2(d) when the strange quark mass m s is set to be zero. 
In the limit x q — > 0, the denominator of the propagator 
is 2fc s • I, which can reach zero when k s and 1 are in the 
same direction. To avoid this singularity, we have taken 
the strange quark mass m s equal to its current mass 90 
MeV in our above numerical calculations. 

There are some uncertainties in the treatment of this 
singularity. To see this, in FIG. 6, we show the results of 
(55 HR as a function of m s for Q 2 = 20 GeV 2 at two £ values 
for two choices of fi. Generally, the results decrease with 
the increase of m s , and in the favored range 70-120 MeV 
for m s , the results may vary within 15%. Lower m s range 
shows greater m s dependence. There is another way to 
avoid the singularity, i.e., to keep the x q suppressed terms 
in the s quark propagator and take x q to be reasonable 
finite instead of zero. When taking x q = 1/6 (approx light 
constitute mass over D meson mass), the numerical results 
reduce nearly by half. 



As indicated by the results above, when measuring 
the strange distribution asymmetry inside a nucleon, one 
should consider the correction 5S^ R (^, Q 2 ) caused by HQR, 
which is comparable to the existent measured value. For 
example, at Q 2 — 20 GeV 2 (about averaged value in 
NuTeV experiment) and for (i = /iq, the HQR correction 
to strangeness asymmetry by integrating (£) over £ 
can be 0.002 approximately. To be comparison, the re- 
cent NLO analysis of the NuTeV dimuon data [17] and 
the global analysis [TI5] present the central value of the 
strangeness asymmetry S~ 0.002. Thus, the HQR could 
not be negligible in the extraction of the strangeness asym- 
metry. 

The HQR correction may enhance the strangeness asym- 
metry by a larger positive value, and large positive strangeness 
asymmetry can help to explain the NuTeV anomaly [6,7, 

The value of the parameter p sm still has some uncer- 
tainty as discussed in Refs. [Mll2"5"] . whereas, the magni- 
tude order of p sm can not be changed and the influence 
from HQR to the measurement of the nucleon strangeness 
asymmetry could not be negligible. Moreover, more accu- 
rate p parameters are needed not only for a better under- 
standing of the HQR effect but also for better determina- 
tion of the strangeness asymmetry inside a nucleon. 



4 Calculation of the LQF-P process 

The diagrams for the LQF-P process in ^-induced CC 
DIS are shown in FIG. 3. The subprocess can contribute 
to the Pa- and F M -induced dimuon cross section difference 
as shown in Eq. (|13p , and thus may also influence the mea- 
surement of the strangeness asymmetry. 
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The calculation of the cross section d&cD for LQF-P 
process can be factorized into the convolution of the LO 
subprocess cross section for light quark q production and 
the fragmentation function (z) of the light quark q into 
a D meson: 



dz da,. 



~-D- 

*fl q^q 



(28) 



where q could be d or s, with D being D + or Df or 
D*+ or D*+. With CP transformation and SU(3) flavor 
symmetry, only two independent light-quark fragmenta- 
tion functions remain: 



D q (z)^Df(z)=Dp(z), 



D* q (z)^D?' + (z) = D^ ' (z). 



(29) 



The fragmentation functions are calculable with pQCD 
30,31,32,33 . For q —> D(cq), the fragmentation function 
can be expressed as |32] 



16tH 



ds 9{s 



{rriq + m c ) 2 m\ 



1 



— ) lim 

Z g°-»oo 



|Mo| ; 



,(30) 



where M and Mo are amplitudes for the D production and 
the LO on shell q production respectively; Let Pd = Pq+Pc 
and k c denote the momenta of the produced D meson and 
the c quark respectively. Then q = Pd + k c is their total 
momentum, and s = q 2 . The variable z is defined by z — 
{P% + Pd)/{q Q + q 3 ) in a frame, where q = (q° , 0, 0, q 3 ). 
In axial gauge, the amplitude corresponding to the dia- 
gram FIG. 3(b) is suppressed and can be neglected [33] . 
Thus only diagram FIG. 3(a) contributes. The formation 
of the bound state can be described by the B-S wave func- 
tions, for the production of 1 S'o state color-singlet D me- 
son, 

R(0) 



XpCSo) 



3^3ttM d 



D) 



(31) 



and for the production of 3 Si state color-singlet D meson, 



X P ( 3 5i) 



i?(0) 



3 % /3^M 



Md + M D ) , 



(32) 



D 



which should appear in the amplitude M of Eq. ([30|) . Here 
i?(0) is the non-relativistic radial wave function at the 
origin for the D meson, and e is the polarization vector of 
the 3 Si state D meson. 

The fragmentation function for light quark into 1 S'o 
state D meson is given by [32 



Dq(z) 



2a s (2m c ) 2 |i?(0)| 2 rz(l - z) 2 

8l7rm 3 [1 - (1 - r)zf 

x [6 - 18(1 - 2r)z + (21 - 74r + 68r 2 )z 2 
-2(1 - r)(6 - 19r + 18r 2 )z 3 
+3(l-r) 2 (l-2r + 2r 2 )z 4 l, (33) 



4.0x10 



3.0x10 



Q 2.0x10 



1.0x10 




Fig. 7. Fragmentation functions D q (z) and D* q {£) for the LQF- 
process. The solid line denotes D q (z), and the dashed line 
enotes D q (z). 



where r = m c jMu- And the fragmentation function for 
light quark into 3 Si state D meson is given by [32] 



D* q (z) 



2a 6 .(2m e ) 2 |fl(0)| 2 rz(l - zf 

27?rm 3 [1 - (1 - r)z} 6 
x [2 - 2(3 - 2r)z + 3(3 - 2r + 4r 2 )z 2 
-2(l-r)(4-r + 2r 2 )z 3 
+ (l-r) 2 (3-2r + 2r 2 )^ 4 ]. (34) 



The value of R(0) can be estimated from the pseu- 
doscalar meson decay constant fp through the relation 



R(0) 



itM d 



h 



(35) 



We take the central values from Ref. [55] for M D + and 
fo+ in calculating R(0), and obtain i?(0) = 0.31 GeV=. 
We take one-loop a s with A = 326 MeV for 4 flavors 
as in CTEQ6L, and obtain a s (2m c ) = 0.255 for m c = 
1.5 GeV. The r value in Eqs. (|33l34|) is evaluated by tak- 
ing M D = 1-87 GeV and M D * = 2.01 GeV[22]. The 
fragmentation functions D q (z) and D*(z) calculated with 
Eqs. (|33|34|) are shown in FIG. 7. If integrating D q (z) and 
D*(z) over z, then D q = / D q {z)dz = 2.01 x 10" 5 and 
D*= J D*(z)dz = 1.77 x 10~ 5 are obtained. 

With the formulas for factorization and the fragmen- 
tation functions Eqs. (|13l28l33l34p . for isoscalar target, the 
'correction' 5S£ CjF _ p from LQF-P to measured strangeness 
asymmetry can be computed: 



5S 



LQF-P 



dxdz- 



7rr 2 x u v (x)+d v (x) 



G 2 F \V cs \ 2 f c B ef 2 
x { ^^^ [D q (z)Br D+ +D q (zyBr D , + }} 
1 



m s \ 2 B e ; {Y ^ BrD++D ^ +)} - (36) 
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Here Y = J x[u v (x) + d v (x)] (1 - yfdx with y = Q 2 /xS < 
1. Y = 0.16 is evaluated from the CTEQ6L parton dis- 
tributions at Q 2 = 20 GcV 2 and E v = 160 GeV, and the 
muonic decay rates are taken to be the central value from 
Ref. [22]. With the parameters given above, the final result 
of c^lqf-p is obtained: 



8S, 



LQF-P 



0.53 x 10" 



(37) 



Such a 'correction' to the strangeness asymmetry from 
the LQF-P process is much smaller than that measured. 
Thus, LQF-P gives little influence in the extraction of the 
nucleon strangeness asymmetry. 

We should note here that as pointed out in Section 
II, LQF-P in neutrino and anti-neutrino DIS may gener- 
ate not only the oppositely charged dimuon events, but 
also the same charged dimuon events and trimuon events 
(not from direct Z?-meson production), thus the fact that 
in neutrino and anti-neutrino DIS experiments either the 
same charged dimuon events or trimuon events are very 
rare is consistant with the small value of 5S^q F _ p as 
shown in Eq. (f3"T|) . 



5 conclusions 

The measurement of the nucleon strangeness asymme- 
try is important for the study of nucleon structure and 
certain related phenomenon. The cross section difference 
between the dimuon production from neutrino and anti- 
neutrino DIS is sensitive observable to the strange distri- 
bution asymmetry. Whereas in this work, we point out two 
types of direct charmed meson production, i.e., HQR and 
LQF-P at order a 2 , which also contribute to oppositely 
charged dimuon production. These processes are not in- 
cluded in the experimental analysis, therefore we further 
study their influence to the measurements of the nucleon 
strangeness asymmetry. With quantitative calculations in 
terms of pQCD, we find that HQR affects the extraction 
of the strange distribution asymmetry with a positive 'cor- 
rection' 5Sft K (£, Q 2 ), and the 'correction' can be so large 
as 5S^ R ~ 10~ 3 . For the other one, LQF-P provides a 
small 'correction' to the measurement, that is of the order 
< 10~ 5 , so that it can be neglected. The influence of HQR 
to the measurement of the nucleon strangeness asymmetry 
from neutrino and anti-neutrino DIS can not be negligi- 
ble, that may provide a positive correction to the present 
value of the strangeness asymmetry, and is also helpful to 
explain the NuTeV anomaly. We think that a reanalysis 
of the strange distribution asymmetry with consideration 
of the direct-D production from HQR process is needed to 
improve the value of the nucleon strangeness asymmetry. 

D-meson directly produced in neutrino and anti-neutrino 
deep inelastic scattering (DIS) contributes to oppositely 
charged dimuon events, which has not been considered in 
experimental analysis so far. Hence we conclude that in 
determining the strangeness asymmetry in a nucleon via 
measuring the dimuon production in neutrino and anti- 
neutrino DIS, the contribution from HQR production of 
the direct D-meson should be deducted precisely. 
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